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Abstract 

We study a supersymmetric phenomenon that can give spectacular signals at the 
LHC: oscillations of neutralinos. Such oscillations can be naturally realised in R- 
symmetric models, where additional fields are introduced as Dirac mass partners of 
gauginos and Higgsinos. Majorana masses for gauginos, necessarily generated from 
anomaly mediation, can create tiny mass splittings between degenerate mass eigen- 
states, causing states produced at the LHC to oscillate between neutralinos and their 
Dirac partner fields. Scenarios where such states decay with displaced vertices can 
lead to striking signatures at the LHC, where the oscillation can be visible directly 
in the distribution of displaced vertex lengths. We elaborate on the theory and LHC 
phenomenology of this feature within a specific scenario of a Higgsino decaying with a 
displaced vertex into a gravitino and a Z boson. 



1 Introduction 



An extensive study of the theoretical and phenomenological aspects of supersymmetry has 
been one of the most prominent facets of high energy phenomenology for several decades, 
and these ideas are now finally being probed by the Large Hadron Collider (LHC). However, 
given the vast richness of theoretical and phenomenological possibilities that supersymmetry 
allows, it is plausible that signals of new physics can show up in spectacular, unexpected 
forms in the coming years. While unexpected signals are exciting for phenomenological 
reasons, they can also be useful in obtaining insights into the details of the underlying 
models. With the LHC rapidly providing new data at energy scales never probed before, the 
study of unexplored phenomenological signatures that could appear at the LHC has become 
extremely important and timely. 

With this motivation, we study a novel phenomenon that can give spectacular signatures 
at the LHC but has remained hitherto unexplored: neutralino oscillations. Oscillations 
within the context of supersymmetry have been discussed in some earlier works; see e.g. 
[H EJ [3] , where the oscillation of a mesino leaves visible signatures in the content of decay 
products at a collider. What we have in mind in this work is a more spectacular signal, 
where the oscillation itself can be directly observed. Such scenarios can be naturally realized 
with Dirac gauginos in the context of R-symmetric models [U El E], where the traditional 
Majorana gaugino mass terms are forbidden by R-symmetry, and new partner fields must be 
introduced to give Dirac masses to the gauginos. However, Majorana masses for gauginos 
will be generated from supergravity effects via anomaly mediation 018]. If these Majorana 
masses are small, they introduce tiny mass splittings between degenerate mass eigenstates, 
enabling a system to oscillate between a gaugino or Higgsino and its partner field. This 
oscillation behavior can leave visible footprints in collider observables; in particular, if the 
states in question are long-lived and decay with displaced vertices within the detector, the 
oscillation can be directly observed in the distribution of lengths of displaced vertices. 

The main purpose of this work is to draw attention to the possibility of this phenomenon 
and explore whether it can be observed at the LHC. For this reason, we will not discuss in 
full detail specific model-building questions or perform full-fledged detector level simulations 
of the signals; such details are premature until hints of such a signal are actually observed. 
We therefore only pursue these directions at a level sufficient to illustrate that the occurrence 
and observation of such a striking phenomenon is feasible at the LHC. 

The paper is organized as follows. We begin by discussing the theoretical motivation 
and details for Dirac gauginos in R-symmetric models in Section [2j followed by a discus- 
sion of the oscillation phenomenon in Section [3j In section [4] we elaborate on the collider 
phenomenology within a specific model of a Higgsino NLSP decaying into a gravitino and 
a Z boson, and discuss the reconstruction of the oscillation feature at the LHC. Section [5] 
contains a broader discussion of other scenarios where oscillations can be realized, and some 
experimental aspects. 
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2 Theoretical Motivation 



We focus on gauginos and Higgsinos in the context of R-symmetric models - in particular, 
the Minimal R-symmetric Supersymmetric Model (MRSSM) |4J - which is well-motivated 
for several reasons. In this section we discuss how Dirac mass terms arise for the gauginos 
and Higgsinos in the MRSSM, and Majorana masses for gauginos are necessarily generated 
via anomaly mediation. We follow the notation in [I]. 

The R-symmetry forbids Majorana mass terms for gauginos, requiring the addition of 
adjoint chiral superfields $j for each gauge group i, with the corresponding fermion carrying 
the opposite R-charge, to accommodate R-invariant Dirac masses. The R-symmetry also 
forbids the \x term, hence also the Higgsino mass, in the Higgs sector, analogously requiring 
the introduction of R u and Rd superfields to partner H u and Hj. 

One can consider the Dirac masses for the gauginos to be elements of a general softly 
broken supersymmetric theory. If SUSY breaking is assumed to originate from hidden sector 
spurions, both F and D-type breaking are allowed, which can be written in terms of the 
spurions as X = 9 2 F and W' a = 6 a D. The W can be thought of as a hidden sector 
U(l)' that acquires a D-term [9]. The Dirac gaugino mass originates from the R-symmetric 
operator involving the D-type spurion [TU] : 



where M is the mediation scale of SUSY breaking from the hidden sector to the visible 
sector, Wi represents the gauge superfield, g is the gaugino, and g' is the corresponding Dirac 
partner. Likewise, in the Higgs sector, the supersymmetric mass term fi u H u R u + fidHdRd 
arises from 



In addition, one also gets couplings of the electroweak $j adjoint chiral superfields to the 
Higgs doublets 



These trilinear terms play a vital role, for instance, in electroweak baryogenesis, and large 
X q (q = u, d) close to the perturbative limit are also favored by a Higgs mass close to 125 
GeV [TT]. These X q couplings affect the mixing in the neutralino sector and are therefore of 
relevance. 

Despite the R-symmetry, Majorana masses for gauginos cannot be completely eliminated. 
It is well known that supersymmetry breaking in a general hidden sector model necessarily 
generates gaugino masses at one- loop as a consequence of the super- Weyl anomaly [7J [8] . This 
anomaly-mediated Majorana mass term, m\, can be obtained through an F-type spurion 




(1) 




(2) 




(3) 



i=B,W 




(4) 
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The consequent Majorana mass for gauginos is given by 



mX = ~2g rm3/2 ' ( ' 

where f3(g 2 ) = dg 2 /din n is the gauge beta function and 1713/2 is the gravitino mass. Recall 
that rrt.3/2 = F/(\/3Mpi); hence the sizes of the anomaly-mediated Majorana masses are 
determined by V~F, the fundamental scale of SUSY breaking. For a bino, for instance, this 
corresponds to [12J 

11a 
mi = 1 2?^^ m 3/2 « 8.9 x 10 6 m 3/2 . 6 

There are no analogous Majorana masses for the R-partners, since the anomaly-mediated 
mass term is given by the running of anomalous dimension and gauge coupling in the field 
strength tensor, which does not exist for the R-partners in the superpotential. Likewise, the 
couplings of the R-partners to Standard Model states are suppressed due to the R-charge 
assignments, and the only fields they couple significantly to are the corresponding R-scalars. 
The R-partners therefore have negligible decay widths of their own when these R-scalars are 
heavier. In the remainder of this paper we assume that this is the case and simply set their 
widths to zero. 

Therefore, the MRSSM has Dirac masses for gauginos, which requires the introduction of 
new partner fields. Majorana masses for gauginos are necessarily introduced from anomaly 
mediation, with their sizes determined by the scale of SUSY breaking. The interplay between 
the Dirac and Majorana masses can have interesting phenomenological consequences. We 
have in mind a Dirac mass of (9(100) GeV and a Majorana mass of O(10~ 5 ) eV. This hierarchy 
opens the possibility of neutralino oscillations, to which we now switch our attention. 



3 Oscillations 

To understand the oscillation phenomenon, consider a gaugino g and its adjoint Dirac partner 
g' . In the (g,g') basis, the Hamiltonian takes the form (see for example [331 OH]) 

<u_( m M - «T/2 m D 
\ m D 

where tum and are the Majorana and Dirac masses for the gauginos respectively, and T 
is the decay width of the gaugino. As discussed in the previous section, we assume that g' 
has a negligible decay width, and no Majorana mass. Note that this negligible decay width 
of g' is also responsible for the vanishing of the decay widths in the off-diagonal entries of 
the Hamiltonian. 

The eigenvalues of this Hamiltonian are 



1*1,2 = 
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For m£)>r, tum, the eigenvalues are approximately 

1 iF ( rriM \ , \ 

W*m D ±-m M --[l±— j . (9) 

The two eigenstates are therefore approximately degenerate with a common mass mo, with 
a small mass splitting tyim between them. In this case, the traditional dimensionless ratios 
x and y are 

Am m M Ar m M 

X = " ^' ^2T = W (10) 
Since we are working in the m£i>r, rriM regime, we set y = from here on. 

Next, let ip\ and ip 2 be the two mass eigenstates, with eigenvalues /^i and fi 2 respectively. 
For infinitesimal mass splitting, one has close to maximal mixing, such that we can write 

^,2^^{~g±~g'), fif,^w-^(Vi±^)- (11) 

First, consider the more general case 

^ a = c al g + c a2 g', (12) 

where a = 1,2. If, at proper time t = 0, the system is ip = x\g + X2fj', its evolution over 
time is given by 

^ (t) = [ e~ iMlt c n (zi c 22 + x 2 c 2 i) - e"^ 2 * c i2 (xj c 21 + x 2 c 22 ) ] g 

+ [ e" iAtl * ci 2 (xi c 22 + x 2 c 2 i) - e"^ 2i c 22 (arj c i2 + x 2 c n ) } g ■ (13) 

Next, consider an interaction that produces a pure gaugino state at t — 0; this is represented 
by (x 1 ,x 2 ) = (1,0). As evident from Eq. (13), at later times this state evolves into a combi- 
nation of the gaugino g and its partner g' . In particular, at time t the gaugino fraction of 
the state is 

(g\ip{t)) = e - * 1 * c n c 22 - e-^ 2i c 12 c 21 w - e - *** 1 * (1 + e~ iniMt ) . (14) 

This shows that the gaugino fraction in the produced state oscillates over time. Since only 
the gaugino, not the Dirac partner, can decay, the probability of decay of the state i/j is also 
correlated to the gaugino fraction, and is proportional to 

l<^)>| 2 = ^e- r *[l + cos(m M t)] (15) 

We are therefore led to a scenario where an oscillation is superimposed on the exponential 
decay traditionally expected of an unstable state. If the state in question travels a measur- 
able distance in the detector before decaying, the oscillation can be observed in spectacular 
fashion: if a large number of such displaced decays are measured, the distribution of the 
displaced vertex lengths will follow Eq. (15), offering a tantalizing collider signature of the 
phenomenon. 
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4 A Specific Scenario: Higgsino Oscillation at the LHC 



We now elaborate on the details of the oscillation behavior and the corresponding signal 
at the LHC within the framework of a specific scenario. We choose a setup that has been 
used in literature to study the reconstruction of displaced vertices in [15]: the decay of a 
Higgsino NLSP into a gravitino LSP and a Z boson. We adopt this process into an MRSSM 
framework together with small anomaly-mediated Majorana masses for the gauginos. The 
Higgsinos have no Majorana masses of their own but mix with the gauginos in the neutralino 
sector and inherit the oscillation behavior discussed in the previous section, while the leptonic 
decay channels of the Z enable clean reconstruction of the positions of the displaced vertices, 
making this an ideal setup in which to study the phenomenon. 



4.1 Parameters 

For simplicity, we assume that the lightest neutralino is a pure up- type Higgsino, H u , which 
can decay into the gravitino LSP and the longitudinal component of the Z-boson. The decay 
length of the Higgsino is given by [151 HB] 



, v . \ 100 GeV \ m 2 7 \ 

ct ^ * r 1 f) x 0.2 mm 16 

l 100TeVy ^ m^o J y m 2 - J 1 ; 

where we have assumed sin 2 /3 ph 1. The decay width depends mainly on the SUSY-breaking 



scale v F and the neutralino mass m^o . 

The mass splitting in this case is slightly involved due to the mixing between the gauginos 
and the Higgsinos in the neutralino sector. Since there is no Majorana mass for Higgsinos, 
the mass splitting between the lightest degenerate eigenstates (combinations of H u and R u ) 
arises through the mixing with the gauginos. The wino mass needs to be above 1 TeV in 
order to satisfy precision electroweak limits on the p parameter [I]. Note that if unification 
of gaugino masses is assumed, this also drives the gluino mass to several TeV, which leads 
to squark- initiated SUSY processes at the LHC as discussed in [17] . The neutralino sector is 
then essentially reduced to the B — H u system (while the wino decouples because it is heavy, 

can be decoupled assuming negligible mixing with this system), for which the neutralino 
mass matrix in the {B, B', H u , R u } basis is 



M: 



N 



( m M(B) m D(B) m z sin W sin (3 \ 

m D(B) 



(17) 



mz sm9w sin/3 — fi u 

V -p u / 

Here we have also ignored the X q couplings from Eq. ^ for simplicity. Nonzero values for 
these couplings do not destroy the oscillation behavior discussed in this section; in particular, 
corrections from their inclusion are negligible for sufficiently heavy R-scalars. These have 
been checked explicitly. 
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gg' 


WW' BB' 


KK 


t*L,RpL,R GG' 


mass (TeV) 


3 


1 0.5 


0.12 


0.6 F/M p i 



Table 1: Mass spectrum for Monte Carlo simulations. 



With the bino sufficiently heavier than the Higgsino (necessary for the lightest neutralino 
to be a pure H u ), the mass splitting between the two lightest neutralino eigenstates is given 
by 

f ^^ Vf^V (18) 

\ D(B) ) \ m D(B) J 

Since R-symmetry is preserved in the absence of Majorana masses, the dependence on the 
bino Majorana mass is understandable. The splitting is additionally suppressed by the 
off-diagonal bino-Higgsino mixing as well as the hierarchy between the bino and Higgsino 
mass scales; this additional suppression is desirable since it translates into larger oscillation 
lengths at the LHC, increasing the likelihood of reconstructing the oscillation feature. Since 
the lightest eigenstate is made up of H u and R u , we also have m^o = // n . 

For the oscillation to be observable at the LHC, both the oscillation scale and the neu- 
tralino decay length in the collider frame must be around mm to m scale, with the decay 
length greater than the oscillation length so that one or more oscillations can be recon- 
structed. From Equations (j6l), (16), and (18), we see that these two scales depend on the 



SUSY breaking scale yF and the neutralino mass m^o. 

We assume the mass spectrum given in Table [l] (the scalars are assumed to be sufficiently 
heavy to be decoupled, and not listed). A mass splitting tiny enough to slow down oscillations 
to the collider scale requires a low scale of SUSY breaking. With 

y/F ~ 3 x 10 2 TeV, (19) 

the mass splitting between the lightest eigenstates is 

Am# u « 10 _5 eV. (20) 

This corresponds to an oscillation scale of cr pa 2ir/Am,jj pa 4 cm. With these scales of 
SUSY breaking and Higgsino mass, the proper lifetime of the neutralino is cr pa 20 cm. 
These decay and oscillation scales are therefore ideal for the observation of this phenomenon 
at the LHC. The above SUSY breaking scale also sets the gravitino mass at the eV scale, 
making it the lightest supersymmetric partner (LSP), as required for the decay channel we 
are interested in. 

All simulations, plots, and discussions in the following sections are specific to these pa- 
rameter values, unless specified otherwise. The approximations and discussions of Equa- 



tions (17), (18) are provided for intuition only; in our runs, we diagonalize the full neutralino 



mass matrix and obtain the eigenvalues numerically. 
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4.2 LHC Phenomenology: Vertex and Oscillation Reconstruction 

We focus on event reconstruction at the ATLAS detector; the general idea is the same at 
CMS. The Higgsino NLSP is promptly produced at the interaction point through the decays 
of stops. Once produced, it travels a macroscopic distance in the collider frame before 
decaying with a displaced vertex into a gravitino and a Z boson. We focus on events where 
the Z boson decays into e + e~, which allows for clean reconstruction of the displaced vertex 
using measurements in the ECAL and tracker. The ECAL can measure energy, timing, and 
pointing information, and this is sufficient to reconstruct the vertex as well as the three- 
momenta of the electron pair. Details of the reconstruction procedure are described in [15] 
and will not be repeated here. The background for such displaced vertex events should be 
negligible, particularly since the lepton pair will also satisfy a Z mass window cut. 

As mentioned earlier, the distribution of displaced vertex lengths should show the oscil- 



lation feature denoted in Eq. (15). There are three main factors that can potentially distort 



the oscillation signal and make it irrevocable; we now discuss each of them in turn. 
4.2.1 Interaction with Detector Material 

The Higgsino travels through material in the detector before decaying. Although hard 
scatterings are unlikely since the Higgsino only interacts weakly, passage through detec- 
tor material can change the oscillation behavior through coherent forward scattering. This 
is analogous to the matter effect or MSW effect for neutrinos [TSJ Q2], and its relevance for 
neutralinos can be estimated by using results from the neutrino MSW effect. 

For neutrinos, interaction with electrons in matter induces an effective potential 

V = y/2G F N e , (21) 

where Gp is the Fermi constant and N e is the electron density. This can be written as [20J 



V = 7.56 • 1(T 14 — ^—r Y e eV (22) 
\g/cm 6 J 

where p is the density of the material and Y e is the number of electrons per nucleon. Since 
neutralinos also interact via the weak force, and the number density of quarks is of the same 



order as N e , Equations (21 ),( 22 ) should provide a reasonable order-of-magnitude estimate 
of the effective potential for the neutralino. With p = 10g/cm 3 and Y e = 1, one obtains 



V ~ 10~ 12 eV. This is several orders of magnitude smaller than the mass splitting in Eq. (20). 
Neutralino interaction with detector material is therefore irrelevant for our study. 

It might appear surprising that the MSW effect is negligible for neutralinos in this case, 
where the mass splitting is only Am m 1CT 5 eV, but can be relevant for neutrinos, where 
the mass splitting is AmJ ~ 10 _5 eV 2 . The main point is that for the MSW effect to be 
negligible, the size of the effective potential should be small compared to the oscillation 
scale, which is given by the difference between the energy eigenvalues. For neutralinos, this 
is given by Am; for neutrinos, it is given by Am^fE, where E is the neutrino energy. This 
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difference comes from the neutralino being massive, which leads to a different subleading 
term in the expansion of the energy eigenvalue. For neutrino energy E rs 10 MeV, we have 
Aml/E rs 10~ 12 eV, which is the same scale as the effective potential V. 

4.2.2 Boost Factors 

The convolution of the neutralino lifetime with the effects of boosting to the reference frame 
of the detector can potentially wash out the oscillation feature. Variations in boost factors 
arise due to convolution with the parton distribution functions (PDFs) for the proton and 
due to kinematics of the various decay processes. The decay probability in the rest frame of 
the Higgsino is 

VQi GZ) oc e- rto [1 + cos (Am t )] . (23) 
In terms of the time measured in the lab frame, this is 

VQi^rGZ) oce- rW7 [l + cos(Amti ab / 7 )], (24) 

and the uncertainty comes from not knowing the value of 7 on an event-by-event basis. 

To overcome this hurdle, we make use of the fact that events with larger 7, corresponding 
to more boosted Higgsinos, are generally correlated with more energetic leptons. Performing 
appropriate cuts on the energies of the leptons should therefore pick out events with approx- 
imately similar boost factors, thereby restoring the oscillation. We verify the efficiency of 
this method with Monte Carlo simulations. To simulate this process, we generate parton 
level events with a Higgsino NLSP decaying into a Z and gravitino in MADGRAPH5 [21] 
for the following process at the 14 TeV LHC: 

pp^t*Ji, ii->th, h^GZ, Z^l + r (25) 

We implement parton-level cuts of \r]\ < 1.5, lepton px > 20 GeV, and an angular separation 
AR > 0.4 between leptons for greater lepton identification efficiency. In order to reduce the 
uncertainty of the angular measurement in the ECAL, we also require each lepton to have 
E > 100 GeV. Note that there might be some combinatoric background when both stop 
decays result in the production of a lepton pair via a Z and the wrong leptons get paired, 
but this should be negligible given the small branching fraction into this final state and the 
Z mass window cut that the lepton pair is required to satisfy. 

The result of performing various cuts on the lepton pair energy is shown in Figure [TJ For 
various cuts, the left column shows the distribution of boost factors and the right column 
shows the corresponding distribution of displaced vertices for the surviving sample (the red 
line corresponds to what would be observed if the process involved a pure decay with no 
oscillations). It is evident that although the convolution with boost factors obfuscates the 
oscillation, a narrower cut on the lepton pair energy, which picks out a narrower distribution 
of 7, can recover the oscillation feature in the distribution of decay lengths. Note that the 
energy cut E > 100 GeV on each lepton also serves as a lower cut on the Higgsino energy, 
and therefore on the value of 7. This serves to eliminate oscillation modes with very short 
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(a) Lepton pair 200 < E < 2000 GeV 



(b) Lepton pair 200 < E < 2000 GeV 




(c) Lepton pair 500 < E < 600 GeV 



0.030 

g 0.025 

° 0.020 

S 0.015 

■8 

-g o.oio 

S-h 

^ 0.005 



0.000 




(e) Lepton pair 550 < E < 560 GeV 
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Figure 1: Distributions of boost factors and the corresponding distributions of displaced 
vertex lengths under different lepton energy cuts. The case of decay without oscillation is 
shown in red. Parton level events are generated using MADGRAPH5 with the cuts and mass 
spectrum as defined in the text. From top to bottom, the fraction of events that survive the 
energy cuts are about 100%, 11% and 1%. 
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Figure 2: Event distribution as a function of f , for lepton pair energies 200 < E < 2000 GeV 
(left) and 500 < E < 2000 GeV (right). 

oscillation lengths that contribute strongly in washing out the oscillation feature, and is 
therefore beneficial. It is also worth noting that the Higgsino lifetime in the lab frame is 
O(ns) (for the 7 factors in Figure [T]), which will not run into problems with the ATLAS 
trigger time delay [22] . 

With a stringent enough cut on the lepton pair energy, the above method is already 
sufficient to recover the oscillation, albeit at the cost of retaining only a small fraction of 
the events. More efficient ways of correcting for the boost factor can be implemented if data 
is scarce. Another approach that uses the correlation between the boost factor and lepton 
pair energy without throwing away events is scaling down the neutralino decay length by 
the lepton pair energy, which approximately corrects for the boost factor. For this purpose, 
define the following quantity to be plotted 

tiit r , . 

f = ~jh> (26) 

where r is the measured length of the displaced vertex, and Ez is the energy of the Z boson as 
reconstructed from the lepton pair energy. Figure |2j shows the improvement with using this 
quantity f instead of the measured length of the displaced vertex, where the improvement 
in reconstructing the oscillation without throwing away a significant number of events is 
evident. 



4.2.3 Uncertainty in Vertex Reconstruction 

There are also uncertainties regarding the reconstruction of the position of the displaced 
vertex, arising from uncertainties in the measured energies and angles at the detector. We 
use the uncertainties as listed in [15] . The largest uncertainty comes from the measurement 
of Odiri the angle of the direction of the lepton relative to the beam as measured in the 
ECAL; it leads to an uncertainty of roughly ~ 2 cm in the position of the displaced vertex. 
We include this uncertainty by smearing the displaced vertex length distribution with a 
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Figure 3: The distribution of displaced vertex lengths with pre-selection and lepton energy 
cut 400 < E < 500 GeV. The distribution before (after) smearing by a Gaussian of width 
a r = 2 cm is shown in orange (blue). 

Gaussian of width 2 cm (the uncertainty in the reconstruction of the primary interaction 
vertex is smaller [15] and can be ignored for this purpose). The distribution of displaced 
vertex distances before and after the smearing are shown in Figure [3} As we can see, the 
effect gets slightly washed out but still persists. The oscillation feature can therefore be 
reconstructed despite measurement uncertainties at the detector. 

4.3 Estimate of data size and parameter space 

We now estimate the amount of data required to observe the oscillation in the scenario 
discussed in this section. For the parameters and mass spectrum specified earlier, and with 
lepton energy cut 400 < E < 500 GeV, we estimate that about 2% of all Z — > events are 
reconstructed and pass all the cuts (this is similar to the efficiency in [15]). To distinguish the 
oscillation feature from the case of a pure exponential decay, we require that a fit with the 
former provide a 3a or better fit to the data than a fit with the latter. For this purpose, we 
only use events with displaced vertex lengths larger than 3cm, which should have negligible 
background contamination from prompt events. With these requirements, we estimate the 
required amount of data to be 

£~ 200fb- 1 . (27) 

This corresponds to ~ 20 displaced vertex events that pass all the cuts. This particular sce- 
nario, at least, is therefore within reach of the LHC within a few years of running. We stress, 
however, that this is an extremely rough estimate and should be interpreted accordingly, and 
can change significantly for a different choice of parameters and cuts. 

Next, we explore the regions of parameter space that can give oscillations observable at 
the LHC. Due to the various experimental constraints discussed in the previous subsection, 
oscillations cannot be visible for all values of the SUSY breaking scale (F) and the higgsino 
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Figure 4: Combinations of SUSY breaking scale, y/F, and higgsino mass, m^, for which the 
oscillation is observable at the LHC with the four constraints discussed in the text, for two 
different bino masses. 

mass {m^). In particular, we place the following requirements: 

• The oscillation wavelength (2n Am- 1 ) must be longer than the precision of the displaced 
vertex measurement to prevent the feature from being smeared out. Here we require 
the wavelength to be larger than 2 cm. 

• The oscillation wavelength must be shorter than 150 cm in order to see at least one 
complete oscillation within the ECAL. 

• The decay length (cr) must be longer than the oscillation wavelength. 

• The decay length must be shorter than ~ 4 x 150 cm, in order to have at least 20% of 
the decays appear in the measurable region. 

Figure [4] shows the parameter space that gives oscillations observable at the LHC based on 
the criteria listed above. We have set 7 = 6 to make these plots. The first two constraints set 
upper and a lower bounds on the SUSY breaking scale y/F, which are reflected in the limits 
in the plots. For a given y/F and bino mass, the other two constraints set upper and lower 
bounds on the Higgsino mass. Since the higgsino production rate is not very sensitive to its 
mass, the observation limit does not depend strongly on m~ h . Improvements in the vertex 
measurement precision or more stringent energy cuts can easily open up more parameter 
space. 

5 Discussions and Conclusions 

Although we have focused exclusively on a single benchmark scenario so far, the oscillation 
phenomenon is more general and can be realized in many other scenarios. We now turn to 
a broader discussion of this phenomenon. 

Broadly speaking, the following criteria must be satisfied for the possibility of oscillations 
at the LHC: 
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1. Two approximately degenerate mass eigenstates, separated by a tiny mass splitting, 
that can be produced in a coherent superposition. 

2. Different interaction behavior (e.g. different decay widths) for the two components that 
the state oscillates between, in order for the oscillation to be observable. 

3. Decay inside the detector for the signal to be detected. 

Any setup that satisfies these requirements can give an interesting oscillation signal at the 
LHC with the right parameters. 

For instance, Dirac gauginos are not necessary for oscillations in the neutralino sector. 
The neutralino sector of the MSSM already contains Majorana masses for the bino and neu- 
tral wino and Dirac masses for the Higgsinos. If the gauginos are several orders of magnitude 
heavier than the Higgsinos, the two lightest neutralino mass eigenstates are combinations 
of H u and H^, degenerate at Higgsino mass /i, separated by a tiny splitting due to differ- 
ent mixing with the gauginos. The mass splitting in this case is 0(n 2 /rrig), which requires 
the gauginos to be at the GUT scale to get oscillations observable at the LHC. A pure H u 
produced in an interaction at the LHC, for instance, will then evolve as a coherent super- 
position of these two mass eigenstates, oscillating between H u and as it travels through 
the detector. If H u and decay into different Standard Model final states, the oscillation 
of the neutralino will manifest itself as variations in the final states produced at different 
distances from the interaction point. 

Likewise, the existence of a gravitino LSP is also not necessary. Neutralinos can de- 
cay with detector-scale displaced vertices in models with R-parity violation (see e.g. [23J, 
where R-parity is broken by bilinear terms hlLH u in the superpotential, and the decays 
are suppressed due to naturally small coefficients and Yukawa couplings, leading to dis- 
placed vertices). If the neutralino sector contains both Dirac and Majorana mass terms, an 
oscillation signal as described in this paper is possible. 

Moreover, oscillations are not confined to the neutralino sector. In the MRSSM, all gaug- 
inos have Dirac masses, and anomaly mediation generates Majorana masses for all gauginos. 
This implies that the charginos and the gluino can also oscillate into the corresponding Dirac 
partners. If the lightest chargino (NLSP) is almost degenerate with the lightest neutralino 
(LSP), the chargino must decay into a three-body final state (e.g. LSP+ virtual W, W into 
two fermions), and can be long lived due to phase space suppression; e.g. [21] has postulated 
charginos with lifetimes ~ 10~ n s. Likewise, the gluino can be long-lived in split-SUSY like 
scenarios where the squarks are heavy and suppress the decay of the gluino. Displaced 
vertices are therefore possible. The problem with observing oscillations with gluinos and 
charginos is that they interact via strong and electromagnetic forces. Long-lived gluinos 
in split-SUSY are expected to get stopped in the detector. Likewise, ATLAS has looked 
for long-lived charginos by searching for disappearing tracks in the tracking volume of the 
detector, where a chargino lighter than 90 GeV with a lifetime between 0.2 and 90 ns has 
been excluded [2H1 [2Z]- Unfortunately, such interactions cause the state to decohere into a 
single mass eigenstate, thereby destroying the oscillation feature. Oscillations with cm-scale 
displaced vertices extending into the body of the detector are therefore possible only with 
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neutralinos, which interact weakly, but not with gluinos or charginos. 

It can nevertheless be possible to uncover evidence of the oscillation even in the absence 
of displaced vertices in some scenarios. Consider a scenario where the decay length is com- 
parable to the oscillation length, so that the state completes only a fraction of an oscillation 
before it decays, but the scale is too small to give rise to displaced vertices, so that the 
decay is prompt. If the two components that the state oscillates between decay to different 
Standard Model final states, the time integrated decay fractions into these different states 
in the presence of an oscillation are different from the fractions that would be expected in 
its absence. Although not as spectacular as the direct observation of oscillations in the dis- 
tribution of displaced vertices, observing this deviation in decay fractions would be a hint of 
an oscillation in effect. 

The observation of oscillations not only provides spectacular signals but can also be 
useful for extracting valuable information about fundamental parameters in the underlying 
theory. In the scenario presented in Section[4| for instance, the oscillation scale is tied to 
the fundamental scale of supersymmetry breaking y~F ; if other parameters of the theory 
have been measured from other processes, the oscillation signal can be used to calculate 
the SUSY breaking scale. The observation of the oscillation in the distribution of displaced 
vertices as discussed in this paper can also be a strong hint of the Dirac nature of gauginos, 
complementing other signatures of Dirac gauginos and R-symmetric models (see e.g. 
[29]). 

Improvements on the experimental front can also improve the chances of observing and 
refining these oscillation signals. The proposed International Linear Collider (ILC) can pro- 
vide significant improvement, in particular by eliminating the uncertainty from proton PDFs, 
which is a significant contribution to the smearing of the signal. Improvements in timing 
techniques at the LHC itself can also significantly improve on measurement uncertainties 
and corrections for boost factors. 

In summary, in this paper we draw attention to the possibility of a neutralino oscillation 
signal that can appear at the LHC. Such a signal can be naturally realized in R-symmetric 
models, where gauginos are Dirac particles but anomaly mediation necessarily generates very 
small Majorana masses for a low SUSY-breaking scale. The oscillation and decay lengths 
can be at the right scales for the oscillation to be observable at the LHC. More detailed 
studies regarding specific model-building issues and detector-level specifics will be required 
if hints of such a feature actually emerge in the data. For the moment, the possibility of 
such a spectacular signal at the LHC, which is reasonably well motivated yet has remained 
hitherto unexplored, is a tantalizing prospect by itself. 
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